Increase of the mean inner Coulomb potential in Au clusters induced by surface 
tension and its implication for electron scattering 
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Electron holography in a transmission electron microscope was applied to measure the phase shift 
Atp induced by Au clusters as a function of the cluster size. Large Atp observed for small Au clusters 
cannot be described by the well-known equation Atp = CsVot {Ce'- interaction constant, Vq: mean 
inner Coulomb potential (MIP) of bulk gold, t: cluster thickness). The rapid increase of the Au 
MIP with decreasing cluster size derived from A<^, can be explained by the compressive strain of 
surface atoms in the cluster. 

PACS numbers: 61.14.Nm, 81.07.Bc, 68.37.Lp 
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Au clusters are considered as prototype material for 
nano-scaled electronic devices and biosensors More- 
over, Au clusters exhibit an exceptional catalytic ac- 
tivity All these potential applications have moti- 
vated numerous studies regarding the properties of Au 
nano-clusters. One fundamental material property is the 
mean inner Coulomb potential (MIP), which plays an 
important role for the quantitative evaluation of experi- 
mental data obtained from electron scattering techniques, 
e.g. transmission electron microscopy (TEM) and elec- 
tron holography (EH). The MIP is the volume- averaged 
electrostatic part of the crystal potential, which can be 
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vifith Planck's constant h, the electron mass and charge 
TO and e, the unit cell volume 17 and the occupation 
number Ui for the atomic species i within the unit cell. 
The important property in Eq.(l) is the atomic scatter- 
ing factor /f'(0) |jy which correlates the MIP with the 
amplitude of the electron wave scattered in forward direc- 
tion. The MIP can be determined by off- axis EH under 
kinematical diffraction conditions according to the rela- 
tion Aif = C^Vbt (Ce-' interaction constant) by measur- 
ing the phase shift Atp between the electron wave passing 
through the sample with a known thickness t and a vac- 
uum reference wave 0/. On the other hand, local TEM 
sample thicknesses can be determined by EH, if precise 
values of the MIP are known, but thus far only MIP val- 
ues for few materials with limited accuracy are available 
0, S [3, 111- For instance, experimental values for 
the MIP of Au between 16.8 and 30.2 V were reported 
0, whereas calculations yield values of 25.0 to 35.9 V 



0, B llii- Moreover, a strong increase of the Au MIP 
up to 45 V was reported for Au clusters deposited on 
Ti02 powder with decreasing cluster size [13]. Recently, 
effective carbon MIP values up to 65 V were reported 
for ultra-thin amorphous carbon (a-C) films compared 
to a bulk value of 9 V [13,] . This indicates that the MIP 
increase for nano-scaled objects could be a general phe- 
nomenon. In this study, we applied EH to determine 1 ) 
the MIP of bulk Au, which corresponds to the MIP of Au 
atoms in the cluster core and 2) the contribution of sur- 
face atoms to the overall MIP of Au clusters to elucidate 
the physical origin of its increase. 

Samples were prepared by low-energy-beam cluster de- 
position of Aun clusters with 10<rt<20 atoms on com- 
mercial a-C substrates, «10 nm thick. Due to the stor- 
age of the sample, a coarsening of the particle sizes oc- 
curs, which leads to Au clusters with diameters D be- 
tween 0.8<Z?<8.0 nm. Off-axis transmission EH was car- 
ried out in a 200 keV Philips CM200 FEG/ST electron 
microscope equipped with a Mollenstedt-Diiker biprism 
in the selected-area aperture holder. Holograms with 
an interference fringe distance of (0.16±0.05) nm and 
a corresponding resolution A£'=(0.32±0.10) nm were 
recorded using a 2048x2048 pixel CCD camera. The 



was reconstructed from the hologram sideband 14] 
by using the DALI program package extended for holo- 
gram reconstruction 15]. For our microscope a constant 
C£;=7.29-10^ rad-(Vm)~^ was determined according to 
Ref. 0. 

Fig. [T^) shows the reconstructed phase shift for an Au 
cluster with D=6.9 nm. The total phase shift of the elec- 
tron wave in the bright region of Fig. [T^) is given by the 
phase shift induced by the cluster and the supporting a-C 
film, whereas the gray background represents the phase 
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shift due to the a-C film only. We use the following proce- 
dure to extract the integrated phase shift induced by the 
Au cluster A(^™^: first, the integration of the phase shift 
along the y-direction of the integration domain (black 
frame in Fig. [T^)) was performed according to: tp™* (x) — 
J^"^ (p{x,y)dy (see Fig. [1]d)). To eliminate the substrate 
contribution fl^i,{x), the background is linearly inter- 
polated (straight line in Fig. Hb)) on both sides of the 
cluster. Finally, At/?™^ is obtained by integration along 
the x-direction: Aip'^l = J^^^^" ((^™*(a;) - ^pHHx)) dx. 
In Fig. [TJ:) Aip'^l of 123 Au clusters is plotted versus the 
radius of cluster projection R. 

To analyze the experimental A(/9™^ we propose a new 
expression for the MIP of Au clusters which distinguishes 
between surface and core atoms. Generally, the atoms in 
a cluster are compressed due to surface tension. Pre- 
viously observed compressive strain in metallic clusters 
(Au, Cu, Ni, Ft) was attributed to the sole contraction 
of the atoms at the cluster surface 1^, which 
is essential for the explanation of the lattice vibrations in 
small particles i20]. Considering that the strain is con- 
fined entirely to the cluster surface, surface atoms are 
under uniform compressive strain, whereas core atoms 
are unstrained. The strain e leads to a decrease of the 
atomic volume for surface atoms to {l+e)^Q,at, where flat 
is the volume of unstrained core atoms, which is identical 
to the volume of bulk atoms. The inverse dependence of 
the MIP on the volume (Vb ~ see Eq.([T])) requires 
accordingly the distinction between the surface and core 
atoms with different atomic volumes. The MIP of Au in 
Au clusters Vj^' can then be expressed by 
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with the number of surface atoms Ns and the total num- 
ber of atoms within the cluster Nt- Vq is the MIP value 
of core atoms, which is equal to the MIP of the bulk mate- 
rial. As expected, Vq^ approaches Vb for bulk material, if 
Ns/Nt^O and e=— 0.02. However, large differences can 
be anticipated for nano-scaled objects, which are charac- 
terised by increased Ns/Nt ratios and significant strain 
e. A theoretical A(/3^^ can then be calculated on the 
basis of Eq. ([2]) as 



t{x,y)^CEVo''n,i, (3) 



where t{x,y) is the cluster thickness at a given position 
within the area of the cluster projection E on the (x,y) 
hologram (image) plane and fid denotes the geometrical 
cluster volume. 

To interpret the measured A(y9™^ using Eq. ([3]) , the clus- 
ter shape and structure needs to be known. According 
to Fig. [T^) the area of the Au-cluster projection corre- 
sponds in a good approximation to disks with radii R 
suggesting that the clusters could be spheres. To verify 
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FIG. 1: a) Gray-scale coded (256 gray levels between black 
(no phase shift) and white (phase shift equal to 2tt rad)) re- 
constructed phase shift obtained for an Au cluster. The frame 
represents the integration domain (see text). The insert shows 
a typical HRTEM image of an Au cluster with M-Dh struc- 
ture; b) O- the experimental phase shift integrated along the 
y-direction. The solid line indicates the linearly interpolated 
background, c) Experimental (Q) and calculated (solid line) 
integrated phase shift of Au clusters Aip''^^{R). The insert 
shows an M-Dh Au cluster with our experimental shape ac- 
cording to [l^ . 



this assumption, high-resolution TEM (HRTEM) images 
of Au clusters with 0.6<-R<3.4 nm were recorded at nor- 
mal illumination (0°) and after tilting the samples by 
±22°. The projected cluster area increases by about 5% 
in the tilted position, which excludes spherical clusters. 
Vertical facets are required with a ratio h/R=OA of the 
vertical facet height (h) with respect to the apparent ra- 
dius R of the cluster projection. 

Calculations indicate that clusters with a pentagonal 
decahedral (Dh) structure variant called Marks-Dh (M- 
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FIG. 2: Uniform compressive strain of surface atoms e{R) for 
M-Dh Au clusters (•: ideal M-Dh clusters). 



Dh) 21] and the face-centered-cubic (fee) truncated oc- 
tahedral (TO) clusters are energetically the most sta- 
ble structures for Aun clusters (50<iV<5000 atoms) 
IgL 13, • A growth study of small Au particles (size 
4-25 nm) produced by gas evaporation in flowing Ar has 
shown that, after condensation of the clusters on a sub- 
strate, the growth of Au particles with icosahcdral (Ih ) 
or Dh structures is favored with respect to fee ones [2J] . 
Accordingly, the Au clusters may have Ih or M-Dh struc- 
tures. But except for the smallest Au clusters, it was 
shown that the Ih structures are energetically noncom- 
petitive as compared with the M-Dh ones [3, [IH] • We 
therefore assume the M-Dh morphology for Au clusters 
in our experimental size range. This assumption is con- 
firmed by HRTEM images showing the typical fivefold 
symmetry (insert in Fig. [T^)). The M-Dh Au clusters 
are characterized by h—OAR, derived from the tilt ex- 
periments, and an apparent diameter of D—2R (see the 
insert in Fig. [TJ;)). The latter relation is assumed to be 
valid because the cluster projections are almost circular. 

To evaluate Eqns.jS]) and (0]), Ns/NriR) and e{R) for 
Au clusters with our shape are required. For the theoreti- 
eal estimation of these properties we distinguish between: 
a) ideal M-Dh Au clusters with closed-shell structures Fj^ / 
and corresponding ideal radii Rid and h) M-Dh Au clus- 
ters without closed-shell structures and intermediate radii 
R. We calculate first Ns/Nt and e for the ideal clusters. 
The total number of Au atoms in ideal M-Dh clusters, 
i.e. the "magic numbers", NT_id can be described by the 
number of atoms n, m and p [lil . 27 1 as shown in the 
insert in Fig. [TJ;). We use only p=2 to approximate the 
observed circular cluster cross section. For a given to, n 
is calculated on the basis of the experimental ratio h/R. 
Then, the difference of NT_id in clusters with (i-hl) and 

(i) closed shells, N^^^ — N^lid ~ -^s_td^ corresponds 
to the number of surface atoms of an ideal cluster with 
(i-hl) shells. 

The following considerations yield estimates for the 
strain e and cluster radii Rid which depend themselves 
on e. Implying that e is confined completely to the 



cluster surface, an ideal M-Dh Au cluster with (i-hl) 

closed shells and radius Rfd^^^ is formed by a) an un- 
strained M-Dh cluster core with (i)-c\osed shells and 
radius Rul consisting of N^^_^j^ core atoms and b) the 
(i-h l)-closed shell, formed by A^gJ^j^^ surface atoms with 
radii Rau_s<Rau- We assume that the unstrained M- 
Dh cluster core consists of unstrained core Au atoms 
with a diameter dAu='2RAu=0 -288 nm, which is given 
by the minimum bond length between Au atoms in 
bulk gold. The strain is proportional to the re- 

duction of the minimum bond lengths between surface 
and core atoms as compared to dAu' e^'^'^^^ = [{Rau + 
Rau.s) — 2i?Ati]/2i?yiu- However, a more useful expres- 
sion for is given by: e{Rf+^^) « [{Rfd^^ - R^) - 
2Rau]/'2'Rau- Rus can be estimated by constructing un- 
strained M-Dh clusters with closed-shell structure us- 
ing only geometrical considerations and Au atoms with 
Rau- To estimate R'f^'^\ we define the atom density of 
clusters as p{Rid)=NT_id/^ci{Rid)- For M-Dh Au clus- 
ters with our experimental shape, an analytical expres- 
sion of flci{Rid) can be derived. The lower boundary 
for Rid (and the upper boundary for e) is calculated 
considering that p{Rid) cannot be larger than the atom 
density of fee bulk Au /96«;fc=58.9 atoms/nm^ by solv- 
ing A^-^*_^]V^cz(-R-d^^'')=P6«ifc- The upper boundary for 
Rid (and the lower boundary for e) is given by pc/=48.0 
atoms/nm^ obtained for a cluster with two closed shells 
from ab initio density functional theory calculations, 
which yields e = —15.9%. The atom density of large 
Au clusters with >10^ atoms is Pd—Pbuik and the strain 
e in these clusters converges towards the strain of sur- 
face atoms in bulk Au between -1 and -2% for low-index 
metal surfaces [l^. The strain e{Rid) in ideal M-Dh Au 
clusters approximated by the average values on the ba- 
sis of pbuik and Pel is plotted in Fig. O The error bars 
represent strain values deduced from phuik and pd- We 
note that Rid differences associated with maximum and 
minimum strain are smaller than the spatial resolution 
in our experiment. 

Values for e and Ns/Nt for M-Dh Au clusters without 
closed shell of surface atoms are obtained by interpola- 
tion between the properties of the ideal ones. The MIP 
value of core Au atoms, identical to the MIP of bulk Au, 
is then calculated by a least-square fit of the experimen- 
tal A(p™^(i?) and calculated phase shift on the basis of 
Eq.(l3]) with Vq as the only free parameter (solid line in 
Fig. [It)). We obtain Vb=(32.2 ± 3.6) V, which agrees 
well with calculated MIP values of bulk Au 0, S] and 
reach, within the error bar, the upper limit of the previ- 
ous experimental MIP data [8|. 

The experimental MIP values of Au in Au clusters 
are estimated from the measured A(p™^(i?) by using 
Vo^'(i?) = Aip'X^iR)/CE^ciiR) (see Fig. E]). The solid 
line in Fig. [3] is calculated using Eq.Q with Vb=32.2 
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FIG. 3: Experimental (Q) and calculated (solid line) Au MIP 
values in M-Dh Au clusters Vq^R). The hatched region rep- 
resents the maximum error of the Vq' values induced by the 
upper and lower strain estimates. 



V, which is the fitted MIP of core Au atoms. Despite 
the scattering of V^'^' for 1.5<R<2.2 nm, Eq.([2]) describes 
adequately the behavior of the MIP. The description is 
particularly reasonable for small clusters with i? « 0.5 
nm where V(f—85 V is obtained. Data points exceed- 
ing the error limits indicated in Fig. [3] can be attributed 
to Au clusters with different h/R ratios, the estimation 
of the background phase shift or a tilted position of the 
cluster on the a-C substrate. The errors associated with 
different possible cluster structures (M-Dh or Ih) on the 
resulting MIP are smaller than the error limits. 

To explain the MIP increase in Au clusters and thin a- 
C films, only effects of changes of the electronic structure, 
permanent electrical charging or adsorbate molecules 
with a lar ge e lectrical dipole have been considered up 
to now [3. 112L Ichikawa et al. attributed the 
increase of the Au MIP to charge transfer from the Au 
clusters to the Ti02 substrate. However, we exclude per- 
sistent charging of the clusters in our study, because the 
Au clusters are deposited on a-C substrate in electrical 
contact with the metallic sample holder. The presence 
of charged and neutral impurities has to be considered. 
The effect of neutral impurities is expected to be negli- 
gible due to the small amount of additional charge den- 
sity. Charged impurities generate an electrostatic poten- 
tial, which contributes to the measured phase shift. How- 
ever, the consisting of our data set indicates that this is 
not a serious problem. We also add that deformations 
of atoms on twin planes in M-Dh clusters cannot be dis- 
tinguished from compressed surface atoms, but this effect 
is estimated to be small compared to surface tension. Fi- 
nally, effects of phonons and surface plasmons in clusters 
on the measured phase shift can be excluded due to the in- 
elastic nature of the scattering process. In our case, the 
strong increase of the Au MIP in Au clusters can be well 



explained by the compressive strain e of surface atoms. 
This effect is present but can not be measured in bulk 
samples, where Ns/Nt is negligible and e is small. In 
contrast, surface contributions dominate and relaxations 
are important in small clusters. 

In summary, we used EH to reconstruct the phase shift 
Aip induced by small Au clusters. Aip is determined by 
the MIP, which can be modeled by distinguishing between 
surface and core atoms with different atomic volumes. 
MIP values of 85 V are measured for Au clusters with 
nm, in contrast to the derived MIP value of Au 
core atoms of Vo=(32.2 ± 3.6) V in good agreement with 
previous calculations for bulk Au. Significant compres- 
sion of surface atoms (reduction of bond length) yields a 
high surface potential contribution to the MIP of clusters. 
An essential additional ingredient of the model, which is 
generally applicable, is the increasing ratio Ns/Nt in 
nano-scaled objects. The effect of the MIP increase in Au 
clusters and in many nano-scaled objects has important 
consequences for the quantification of electron scattering 
data in general because the amplitude of the electron wave 
scattered in forward direction can be much higher than 
expected from atomic scattering factors given e.g. in f^]. 
Modified scattering factors should be considered for the 
evaluation of images and diffraction patterns of nano- 
scaled objects from techniques based on electron scatter- 
ing. 
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